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Abstract. 

We assess the reliability of the one-crossing approximation (OCA) approach in 
quantitative description of the Mott transition in the framework of the dynamical 
mean field theory (DMFT). The OCA approach has been applied in the conjunction 
with DMFT to a number of heavy-fermion, actinide, transition metal compounds, 
and nanoscale systems. However, several recent studies in the framework of impurity 
models pointed out to serious deficiencies of OCA and raised questions regarding its 
reliability. Here we consider a single band Hubbard model on the Bethe lattice at finite 
temperatures and compare the results of OCA to those of a numerically exact quantum 
Monte Carlo (QMC) method. The temperature-local repulsion U phase diagram for 
the particle-hole symmetric case obtained by OCA is in good agreement with that of 
QMC, with the metal-insulator transition captured very well. We find, however, that 
the insulator to metal transition is shifted to higher values of U and, simultaneously, 
correlations in the metallic phase are significantly overestimated. This counter-intuitive 
behavior is due to simultaneous underestimations of the Kondo scale in the metallic 
phase and the size of the insulating gap. We trace the underestimation of the insulating 
gap to that of the second moment of the high-frequency expansion of the impurity 
spectral density. Calculations for the system away from the particle-hole symmetric 
case are also presented and discussed. 


PACS numbers: 73.23.-b, 71.10.Hf, 75.20.Hr 
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1. Introduction 

In the past years, many efforts have been devoted to the implementation of calculation 
techniques to describe the electronic structure of strongly correlated complex materials. 
This is a complicated and challenging task in view of the many degrees of freedom 
involved. One of the most successful approaches in this direction was the implementation 
of the dynamical-mean held theory (DMFT) [1, 2, 3]. Numerically, the most challenging 
part of DMFT is the solution of the Anderson impurity model [4] within the DMFT 
self-consistent loop that maps the lattice problem into a single impurity one. 

There are two well-known numerically exact techniques to solve this impurity model, 
namely, the quantum Monte Carlo (QMC) in its Hirsch-Fye (HF-QMC) or continuous 
time (CT-QMC) versions [5, 6], and the numerical renormalization group (NRG) [7, 8]. 
Recently, a substantial technical progress [9] has been achieved in both approaches. On 
one hand, the advent of continuous-time quantum Monte Carlo methods [10] eliminated 
the time discretization error, inherent to the HF-QMC, and extended the range of 
applicability of QMC to much lower temperatures and realistic Coulomb repulsion 
vertices. On the other hand, very fast implementations of NRG applied to multi-band 
systems has been developped using Abelian and non-Abelian symmetries on a generic 
level [11]. 

In spite of recent technical improvements, those exact methods still encounter 
certain difficulties. QMC solvers suffer from the well known ’fermion sign problem’, 
which can be especially severe when the degeneracy of the correlated shell is large and 
signihcant off-diagonal terms are present in the hybridization function. Moreover, QMC 
calculations are carried out in the imaginary-time domain and an analytic continuation 
is required to obtain real-energy spectral functions from QMC data. The NRG approach 
becomes computationally expensive in multiorbital cases with broken orbital symmetries 
(for instance, when interactions, like pair-hopping, prohibit the use of symmetries that 
reduce the size of the matrix to be diagonalized [12], leading to an exponential increase 
of the Hilbert space). Because of these limitations the necessity to have faster and 
reliable impurity solvers is evident. 

Hence, several approximate schemes have been proposed for solving the 
DMFT impurity problem, like the local moment approximation (LMA)[13], iterative 
perturbation theory (IPT)[14], exact diagonalization[15], rotationally invariant slave 
bosons [16], conserving diagrammatic approximations based on self-consistent 
hybridization expansion (SCH) [17], among others. 

Regarding the SCH, the non-crossing approximation (NCA) [18] represents the 
simplest family of these self consistent treatments and provides an accurate calculation 
of the impurity Green function, as well as many other properties, when the Coulomb 
repulsion is taken to be large enough as compared with the other energy scales involved 
in the problem. However, when more than one charge fluctuation needs to be included 
{N N — 1 and A^ ^ A^ -f 1, being N the impurity valence), NCA has failed to give 
the correct Kondo scale (Tx). The next leading order in the self consistent expansion. 
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that partially solves this pathology, is often known as the one-crossing approximation, 
OCA [19, 20, 21]. Within this extended formalism other classes of problems have been 
investigated [22, 23, 24], Among them, its major application is in the context of the 
dynamical mean held theory as an impurity solver [3]. 

In particular, the OCA solver has the advantage of being formulated at the real 
frequency axis and it gives the correct order of magnitude for the Kondo scale of 
the impurity problem. It successfully captures the correct temperature dependence 
of transport properties of a single impurity level [21], and it has been employed as the 
DMFT impurity solver in a search for signatures of a non-Fermi liquid behavior in the 
Hubbard model with van Hove singularities [24]. Furthermore, it has been generalized 
to an arbitrary number of orbitals and interactions [23]. Multiorbital generalization of 
OCA were employed in a study of the itinerant and local-moment magnetism in the 
three-band Hubbard model [27]. In combination with ab-initio-\-DMFT calculations, 
the OCA solver has been applied to real strongly correlated materials, for example, to 
heavy-fermion compounds [3, 23, 25, 22]. 

However, the OCA solver has also several limitations. It cannot be applied to 
arbitrary low temperatures due to violations of the Fermi-liquid properties (in the 
impurity model, OCA works well for T > O.ITk) [21, 28, 29, 30], and it also violates 
the sum rules for the coefficients of the high-frequency expansion of the self-energy 
[31]. While the former pathology can be controlled by restricting its application to 
high enough temperatures, the later one is intrinsic and will always be present. As has 
been pointed out recently, the OCA method is more accurate in the strongly-correlated 
limit [31], and it describes the insulating phase particularly well [32]. It has also been 
shown that OCA overestimates the correlations in th metallic phase and it has been 
conjectured that this overestimation of correlation effects reffects the fact that the OCA 
tends to favor the insulating state. 

One important issue that has not been studied up to date is the actual quantitative 
performance of the OCA solver within DMFT in describing the metal-insulator Mott 
transition [33]. Hence, we address this issue in the present work by calculating the critical 
Uc values for the Mott transition within DMFT as a function of temperature using OCA 
as the impurity solver, and comparing them with the corresponding ones obtained with 
the CT-QMC. We have also compared the DMFT local self-energies obtained within 
the two approaches as well as the corresponding quasi-particle effective masses in the 
metallic phase. Our calculations have been carried out for the single band Hubbard 
model with a semicircular non-interacting density of states. 

Our main conclusion is that the OCA metal-to-insulator transition for the particle- 
hole symmetric case is in remarkably good agreement with that of CT-QMC. However, 
we find that insulator-to-metal transition is shifted to higher values of U despite the 
fact the correlations of the metallic phase are overestimated. This counter-intuitive 
behaviour is explained as a combination of two factors: the underestimation of the 
effective Kondo temperature in the metallic phase and the underestimation of the gap 
in the insulating one. The fact that OCA underestimates the gap in the insulating 
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regime comes out from an analysis of the high-frequncy expansion sum rules of the 
Green function. Our results are in contradiction to the conjecture of OCA favoring the 
insulating phase. We show that although OCA overestimates the strength of correlations 
in the metallic phase, it does not favor the insulating one because the critical values of 
the metal-to-insulator transition are very well captured. 

We have also study the same model in the non-symmetric case, obtaining similar 
agreement between both techniques. We verify that the OCA approximation does not 
violate the Friedel sum rule in the metallic phase for the range of temperatnres of the 
obtained phase diagram, and that the interacting part of the OCA self-energy always 
remains causal. 

The paper is organized as follows: we describe the theoretical formalism in section 
2, we present the numerical results for the particle-hole symmetric case in section 3.1, 
we discuss the results obtained for the system away from half-hlling in section 3.2 and 
hnally we conclude in section 4. 

2. Model and Formalism 

We start with the single-band Hubbard Hamiltonian, 

H = ^ ( 1 ) 

^ <ij>cr i 

where the hrst term is the kinetic energy, t is the hopping between nearest neighbors 
on a lattice, 2 ; is the coordination number, and U is the energy of the on-site Coulomb 
repulsion. The operator creates an electron with the spin a on the site i and Ujo- = 
cl^Cia- We use the semicircular non-interacting density of states N{u) = 
leu I < 2t corresponding to a Bethe lattice with coordination z —)■ 00 , for which the 
DMFT approximation becomes exact. In the following we use the half bandwidth as 
our unit of energy D = 2t = 1. 

We solve the Hamiltonian [1] by means of DMFT, which maps the lattice model 
onto a single-impurity Anderson one within a self-consistent cycle. The hopping between 
the impurity and the conduction band, I 4 , dehnes the hybridization function for the 
single-impurity problem F(ia;) = /{ioj — Ck), where Ck is the conduction energy of 

the impurity model. Within the DMFT and in the case of the Bethe lattice, the DMFT 
hybridization function is given by the self-consistency condition F(ia;) = t‘^G[r{iu)], 
where G{u) is the local Green function obtained from the impnrity model. 

Starting from the metallic non-interacting solution of the model, the system turns 
into an insulator for large enough values of the Coulomb repulsion U due to the vanishing 
of the quasiparticle weight. The value of U = Uc 2 dehnes this transition. On the other 
hand, starting from an insulating solution, the systems turns metallic due to the collapse 
of the gap between the Hubbard bands, for U < Ud, with Ud < Uc 2 when T is lower 
than the second-order end point of the hrst-order Mott transition T^. The critical values 
Ud <U < Uc 2 as function of the temperature T determine a phase diagram. 
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The phase diagram of the Mott transition for the present model have been 
previously obtained using the QMC [34, 35, 36], IPT [1], exact diagonalization 
[37, 38, 39], and NRG [8] impurity solvers. The determination of the exact boundaries of 
the coexistence region has previously required a signihcant effort due to their sensitivity 
to calculational parameters, as well as due to the critical slowing down of the DMFT 
convergence close to those boundaries [35]. Hence, we have employed up to 220 DMFT 
cycles for each point in the {U,T} space and used a dense mesh along the U axis with the 
spacing between U values down to 0.005 in the vicinity of the Ud line. We have used 
the CT-QMC implementation provided by the TRIQS package [40, 41]. The DMFT 
impurity problem has been solved by CT-QMC using 10® CT-QMC moves with each 
200 moves followed by a measurement. The resulting CT-QMC phase diagram is in 
agreement with the extensive HF-QMC calculations of Bliimer [36]. Within the OCA 
solver we have used the procedure described by Hettler et al. for regularizing the spectral 
functions [42] and the numerical convolution sketched in Ref. [22] when computing the 
self-energies and the Green function. 

3. Numerical Results 

In this section, we present the numerical results obtained using the OCA solver for the 
DMFT loop and a detailed comparison with CT-QMC calculations. 

3.1. Mott transition for the particle-hole symmetric case 

In order to get the critical values Uci{T) and Uc 2 {T) for a given temperature T within 
the OCA solver, we take advantage of its self-consistent nature building an external 
loop running in the U values. Starting from a metallic solution we slowly increase U by 
6U retaining the previous ionic self-energies and Green function as the initial guess for 
the following U -|- 6U DMFT cycle, until an insulator solution is reached, and then we 
decrease U by —6U until we go back to the initial U. 

In Fig. [la] we show the spectral weight at the zero-frequency, A{u) = 0) = 
—Zm[G{u = 0)]/7r, as a function of U for an inverse temperature f3 = 80. We show both 
the increasing U results from the metallic to insulator solutions as well as the decreasing 
ones. An hysteresis curve is formed, giving rise to two different critical values, Uci{T) 
and Uc 2 {T). We dehne these critical values following the criteria given in Ref. [8], from 
the t/-value for which \A'{u: = 0)| reaches its maximum intensity. 

In Fig. [lb] we show the variation of the quasiparticle weights, Z = [1 — 

as a function of f/ for = 80. In order to compared with CT-QMC, 
we hrst obtain the interacting part of the OCA self-energy D(a;), removing the non¬ 
interacting offset given by the hybridization term. Secondly, from a Hilbert transform 
of we compute the corresponding self-energy in the Matsubara domain, 

1 f , ImZico) 

— I du -. 

IUJ„ — OJ 


Z{iuJn) 


71 


( 2 ) 
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Figure 1. (Color online) a). Spectral weight A{ui = 0) for the inverse temperature 
/3 = 80 as a function of U both for increasing (black lines, squares) and decreasing 
(red lines, circles) U values. The CT-QMC (OCA) data are displayed with the solid 
(dashed) lines and empty (filled) symbols, respectively, b). The quasi-particle residue 
Z as function of U for the same temperature. The notation is the same as in panel a). 

Finally, we approximate the derivative U=o = by a cubic 

fitting of the hrst four Matsubara’s frequencies of ImT,{iuJn). 

Although the vanishing of Z dehnes the critical value Uc only at zero temperature 
[8] , it has been used as a common criteria even for hnite temperatures (see for instance 
Ref. [43]) From Fig. [1] it can be seen that both approaches (from A{u: = 0) or from 
Z) dehne the same energy scales for Uci and Uc 2 - More importantly, the OCA critical 
17-values are in a reasonable agreement with the CT-QMC ones. While the OCA value 
for Uc 2 is obtained within an error of less than 0.5% with respect to the CT-QMC one, 
the calculated Ud is larger than the CT-QMC one by around 3% . We will discuss the 
origin of this discrepancy for Ud later in this section. 

It is important to remark that the OCA values of Z in the metallic region, i.e. 
U < Ud,Uc 2 , are smaller than the CT-QMC ones. The same behavior was found by 
Schmitt et al. [27] using OCA for a body-centered-cubic lattice in comparison with NRG 
calculations. While OCA gives the correct low energy scale for the impurity model, this 
energy scale is still slightly underestimated [20], and therefore within OCA the system 
feels a larger effective Coulomb repulsion giving rise to a reduced quasiparticle weight. 
However, it is important to remark that the underestimation of Z is less important close 
to the transition. 

In Fig. (2) we show the imaginary part of the self-energy in the imaginary frequency 
domain for the increasing U regime at /3 = 60 and for two different values of U, one 
below and one above Uc 2 , U = 2.3 and U = 2.4. As it can be observed from this plot, 
for the metallic case, OCA overestimates the absolute magnitude of the self-energy at 
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low frequencies. Similarly to the underestimation of the quasiparticle weigth at low 
temperatures described above, this behavior of can be also understood as 

arising due to an effectively larger value of U. On the other hand, in the insulating 
region the agreement between OCA and CT-QMC is remarkable. We found that for a 
correct comparison between the two techniques it was very important to have the same 
degree of precision of the convergence criterion of the DMFT loops, especially for points 
close to the Mott transtion. For large frequencies, an additional test can be done using 
the sum rules that T,{iojn) should satisfy. 



Figure 2. (Color online) Comparison of the imaginary part of the self-energy as a 
function of the Matsubara frequency between OCA and CT-QMC at /? = 60 for two 
different values of U, one below the Uc 2 and the other one above. The inset shows 
the imaginary part of the OCA self-energy scaled by a;„. The dashed and solid lines 
indicate their expected theoretical values given by the high frequency expansion sum 
rule, Si = —1/^/4. 

In the inset of Fig. (2) we plot the imaginary part of the OCA self-energy scaled by 
Un for U = 2.3 < Uc 2 and U = 2.4 > Uc 2 , together with the exact coefficient Si for each 
U, that corresponds to the hrst moment in the self-energy high frequency expansion. 
Si = / ^ Im'Z{uj), and that determines the asymptotic l/a;„ behavior. In Ref.[31], 
Riiegg et al. have calculated the exact value expected for Si, being Si = —17^/4 for the 
symmetric case [44]. For the parameters shown in Fig. (2), we obtain a deviation of the 
OCA Si coefficient of the order of 5% in the metallic phase, while in the insulator one 
the error is reduced to less than 2%. 

In what follows we discuss the phase diagram of the Mott transition. In Fig. [3] 
we show the T vs. U diagram with the calculated Uc^ and Uc 2 obtained from the 
zero-frequency spectral function A{uj = 0) (upper panel), as well as the quasiparticle 
residue Z (lower panel). The general trend of the critical Uc{T) obtained by OCA is 
in reasonable agreement with the corresponding CT-QMC ones. Even though a very 
well dehned coexistence region is captured by OCA, this coexistence region is reduced 
with respect to the CT-QMC one. While the agreement is remarkable for the Uc 2 {T) 











Reliability of the one-crossing approximation in describing the Mott transition 
transition, the Uci{T) valnes are slightly shifted to higher energies in OCA. 



Figure 3. (Color online) The T vs. U phase diagram of the Mott transition obtained 
from the zero-frequency spectral function A(uj = 0) (upper panel) and the quasiparticle 
residue Z (lower panel). The inset in the upper panel indicates the phase diagram 
obtained using the finite-U NCA in the symmetric case as the impurity solver. 

Regarding the critical temperatnre (T^) below which two different spinodal lines 
dehne the coexistence region of the insulating and metalic regimes of the Mott transition, 
OCA gives Tc ~ 0.02 in reasonable agreement with the CT-QMC T^. ~ 0.025. The 
slight underestimation of is a consequence of the corresponding underestimation of 
Tk by OCA at the effective impurity level. For comparison, we also include in the 
inset of the upper panel of Fig. [3] the hnite-U NCA phase diagram for the particle-hole 
symmetric case. We stress here that this simple approximation severely underestimates 
all the energy scales involved, Tc and both Uci(T) and Uc 2 (T), as a consequence of the 
underestimated Kondo scale. On the other hand, we want to mention here that the IPT 
results [1, 8] are considerably shifted to higher energies overestimating both, Uci{T) and 
Uc 2 {T), due to the exaggerated overestimation of the Kondo scale at the impurity level. 

Despite its approximate nature, the coexistence region given by OCA is in the 
correct energy range and the critical temperature Tc is in very good agreement with the 
CT-QMC results. We want to remark that for the whole range of temperatures studied 
in the presented phase diagram, the OCA self-energy remains causal, that is, ImT,{iUn) 
is negative. For very low temperatures (T < 1/500 ~ 0.1 it can turn positive 

signaling the breakdown of the approximation. 

We turn now to the discussion regarding the slight overestimation of Ud that can 
be observed in Fig. (3). While the value of Uc 2 is given by the critical U for which the 
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quasiparticle weight at zero frequency vanishes, the Ud is related to the corresponding 
U for which the Hubbard bands collapse and the gap in the spectral function is closed. 
We found that the size of the gap in the insulator regime given by OCA is somewhat 
underestimated and therefore it closes for a larger value of U than expected for CT- 
QMC. This statement follows from an analysis of the high frequency expansion of the 
local Green function. As described in Ref. [31], the high frequency expansion in the 
imaginary domain of G{iUn) is given by 


G{iUJn) 


oo 


E 


Mk-i 


( 3 ) 


where, in the spectral representation of the Green function, the coefficients are related 
to the moments of the spectral density as doo u^A^u) ]; . Exact relations 

for the coefficients can be found from thermodynamic expectation values [31]: Mq = 1, 
Ml = ed + Und/2 (0 at half hlling), and M 2 = Ai-\-U{2ed-\-U)nd/2. Here, and 

are the energy level and total occupancy of the effective Anderson model. Mq and Mi 
are related to the normalization and parity of A{u}) so that they are exactly reproduced 
by OCA. 

Regarding the coefficient M 2 , the parameter Aq represents the zero moment in the 
hybridization high frequency expansion, Aq = —^J^du XmT{u)) = ^ du A(ci;), 
where A(a;) = 7rV‘^pc{uj), and pc is the conduction density of states. Using the self- 
consistency condition r(za;) = t^G\r{iuj)] for the present case of the Bethe lattice, 
we arrive to the following relation: A(a;) = nt^A^uj) = Therefore, Aq = 


d2 


duj A{uj) = Taking into account that for the symmetric situation 2ed-\-U = 0 


and Ml = 0, the coefficient M 2 reads 
f/2 


( 4 ) 


The second moment M 2 of the spectral function contains indirect information about 
the size of the Mott gap. In fact, it carries information about the center position and 
width of each Hubbard band. For instance, in the simplest case in which the Hubbard 
bands have the semicircular shape centered at iwo and width D, the second moment 
becomes M 2 = Uq D‘^/4 by comparing with Eq.(4), one can infer that uq = U/2. In 
this simple picture, the gap is opened when U is larger than 2D and the size of the 
gap is of the order oi 5 = U — 2D. In Fig. (4), we show the spectral density in the 
insulating region when decreasing the Goulomb repulsion from U = 3 to U = 2.6. It 
can be observed that the gap is continuosly closed when U is lowered until the critical 
value Uci is reached. In the inset of Fig. (4), we show the values of doo co^A(lo) 

(squares), which represents the ratio of the second moment obtained within OGA and its 
exact value from Eq.(4), as a function of U and its deviation from the unity (solid line). 
It can be seen that OGA underestimates the second moment of the spectral function by 
~ 15%. 


f With our notation the moments Mk are equal to the coefficients Ck+i defined in Ref. [31] 
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Figure 4. (Color online) Spectral density in the insulating region when decreasing 
the Coulomb repulsion from U = 3 to U = 2.6. The inset shows the ratio of the second 
moment obtained within OCA and its exact value from Eq.(4) (squares) as a function 
of U and its deviation from the unity (solid line). 


Unfortunately, the center position and width of each Hubbard band enter in 
a combination within M 2 and we cannot know from this coefficient alone, if OCA 
underestimated the center position or width or even both. However, an underestimation 
in both quantities bring about a reduccion of the gap that gives rise to larger values of 
Uci as compared with the exact CT-QMC ones. 

3.2. Non-symmetric case 

In this subsection, we compare the calculations done by OCA and CT-QMC for the 
one band Hubbard model on the Bethe lattice away from half-hlling. We consider 
2.5 < U < 5.0 and the impurity level of the effective Anderson model at = — ^ -|- A/i, 
with A/i=-1.0 and (3 = 60. 

In Fig. 5, the spectral densities calculated by OCA for different values of U are 
shown. One may see that for the smallest value of (7, the system is metallic with a 
large quasiparticle resonance that overlaps with the upper Hubbard band giving rise to 
large charge fluctuations pertaining to a mixed valence regime. In the other extreme, 
for the largest value of U, the systems is an insulator with the Hubbard bands located 
symmetrically with respect to A/r. The value of the gap in this case is of the order of 
2D. In order to be able to describe accurately solutions with large gaps, we implemented 
a three-centered logarithmic mesh. 

By integrating A{u) weigthed by the Fermi function for the corresponding 
temperature, we obtained the local occupancies in a very good agreement with the 
CT-QMC ones. It is not obvious that this quantity can be correctly evaluated within 
approximate analytical solvers. Hence, the fact that it is captured within OCA is 
important for the applicability of the method to non-symmetric cases. 
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In the inset of Fig. 5 we show A{u = 0) as a fnnction of U in comparison with 
CT-QMC. One sees that both the OCA and CT-QMC indicate that the system turns 
an insulator for U > 4.5. For this level of doping there is no coexistence region and the 
OCA critical value Uc agrees with the CT-QMC one within a 5%. 



Figure 5. (Color online) Spectral density A{uj) calculated by OCA for a non- 
symmetric case taking 2.5 < U < 5.0 and an energy shift of-1.0 from the corresponding 
symmetric case for each value of U. The inverse temperature is /3 = 60. In the inset 
we show A(uj = 0) as a function of U. The CT-QMC (OCA) data are displayed with 
the solid (dashed) lines and empty (filled) symbols, respectively. 

Overall, we show that OCA also gives a very reasonable description of the Mott 
metal-insulating transition for the Hubbard model away from half-hlling. 


4. Summary and conclusions 

The self-consistent hybridization expansions in their different forms (NCA, OCA, 
symmetric £nite-U NCA, etc) have been widely used not only in the context of the 
impurity problem, but also in the framework of DMFT applied both to different lattice 
models and realistic cases, describing strongly correlated materials from hrst-principles. 
However, to the best of our knowledge, a detailed and quantitative study of the Mott 
transition, one of the essential problems of strongly correlated systems, has not been 
carried out up to now with these kind of approximate techniques. 

In this work, we asses the reliability of OCA impurity solver in the context of the 
DMFT method to describe the Mott metal-insulator for the one band Hubbard model 
in the Bethe lattice at half-hlling within DMFT. We present the temperature-local 
repulsion U phase diagram in comparison with the numerically exact CT-QMC. We 
show that OCA can provide a very good quantitative description of the metal-insulator 
transition of the present model. We obtain the metal-to-insulator transition, Uc^, within 
an error of less than 0.5% while the insulator-to-metal Uc^ values are shifted to higher U 
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(about a 3%) with respect to the CT-QMC one. We explain the overestimation of Uc^ 
from an analysis of the second moment of the spectral density, M 2 . We hnd that the 
expected theoretical value for M 2 is underestimated by OCA. Since M 2 is equal to the 
second moment of the spectral weight, we infer that the size of the gap in the insulating 
phase is also underestimated so that the Hubbard bands collapse for higher values of U 
than for CT-QMC. 

Aside from the Mott transition itself, we conhrm previous results[31, 32] regarding 
the better performance of OCA in the insulating phase than in the metallic one. The 
high-frequency sum rules for the imaginary part of T,{iu) are obtained reasonably well 
in both phases, with the deviation in the insulating case being a bit smaller than in 
the metallic one. On the other hand, in the small frequency region the correlations 
are overestimated in the metallic case. This effect is also apparent in the value of 
the quasiparticle weigth that is underestimated by OCA, specially far away from the 
transition. This overestimation of the correlations in the metallic phase does not imply 
that OCA favors the insulating state, as has been previously stated in Ref. [32], since 
we show the transition U is well reproduced, especially the Uc^ values. Furthermore, we 
show that the gap of the insulating phase is underestimated by OCA. 

Finally, we study the perfomance of OCA for a non-symmetric case obtaining an 
overall reasonable agreement with CT-QMC, and a very similar critical value of U for 
the Mott transition at the considered temperature. The study of non-symmetric cases 
are particularly relevant for applications to real materials. 

Despite the above mentioned deviations of OCA from exacts results, we are not 
aware of any other approximated technique yielding a phase-diagram with this level of 
agreement with numerically-exact many-body methods. 
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